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ABSTRACT 







The propagation of a simple acoustic transient in an isoveLocity 
water layer was investigated for the purpose of studying the region 
of validity of a suggested correlation between exact solutions in 
ducts and solutions for propagation in a layer. The Laplace trans- 
form method was used to obtain approximate solutions for the acoustic 
pressure, particle velocity, and particle displacement resulting 
from a step- function input in velocity. The radiation impedance 
characteristics of a Mylar transducer were investigated to study 
their effect on the received waveforms. Computer programs were 
written to evaluate and graph the resulting waveforms. The wave- 
forms resulting from a step- function input in velocity were observed 
and compared with the predicted waveforms. With the use of a family 
of Mylar transducers, a somewhat doubtful correlation was obtained 
in the region of validity of the theoretical solutions. Possible 
causes of the disagreement between the theory and the experiment 
were investigated. It was found that the theory is believed to be 
correct, that the experiment was not accurate enough, and that 
further investigations should be conducted to study the assumed 
pressure-release characteristics of \ inch polyethylene. 
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1. INTRODUCTION 



Although most of the current investigations of the propagation 
of energy in waveguides have been done for electromagnetic cases . 
some have considered the propagation of pulsed acoustic energy in 
waveguides^ The dispersion of an ultrasonic pulse propagated 

through a waveguide of circular cross-section filled with water 
has been used for the approximate determination of the frequency 
spectrum of the original pulse^\ J a M. Proud et al /"^ have shown 
that the dispersion of a rectangular pulse along an acoustic wave- 
guide of rectangular cross-section can be used to calibrate the wave- 
guide so that it can be used to predict the effects of similar 
dispersion on the transmission of pulses of other shapes that are 
not amenable to actual computation « It has been shown, by considering 
pulse compression in an acoustic waveguide, that a frequency modulated 
input pulse can be compressed into an output pulse of shorter duration 

( 3 ) 

and of higher peak amplitude a In addition, studies of transient 
propagation in layered media have been made v ’ / 0 Many of these 
studies have used approximation techniques, such as the method of 
stationary phase, to solve the transient propagation problem. J. R. 

(g) 

Wait has presented a summary of these techniques useful for cal- 
culating the resultant propagated waveform. Using the Laplace trans- 

(5 9 10 ) 

form technique, others 9 y have investigated transient radiation 
and scattering from cylindrical bodies. 

Recently, exact solutions for the propagation of simple transients 
in waveguides have been published^’ ^ ^ . It has been suggested 
that there is a correlation within a certain region of validity between 
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exact solutions in ducts and solutions for the propagation of 

. . n , (5,14) 

acoustic transients in isovelocity layers 

The purpose of this research is to extend investigation of 
the region of validity of the suggested correlation and of the 
restrictions under which the solutions describing the transient 
response of the layer are valid. 
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2. THEORETICAL DEVELOPMENT 



Development of Solutions 

The response of an isovelocity layer to simple transient 
excitations will be developed using the following assumptions: 

(a) cylindrical velocity- source of radius a, 

(b) constant layer depth d, and 

(c) pressure-release upper and lower boundaries. 

The wave equation for acoustic velocity-potential in cylindrical 
coordinates, under the assumption of radially- symmetric motion, is 

( & * + i ♦ & - k fe) § -- o 2-1 

where r is the radial direction and z is the transverse direction. 



The convenient separation of variables 



casts Eq. 2.1 into 



[ 






t * 



1. A 
r 






[& +(^) z ] if*)- o 



2.2 



where the radial component of the particle velocity becomes 

Ur - Z(i) 

and the acoustic pressure becomes 

For pressure-release top and bottom, we have p = 0 at z = 0 



and z = d; therefore, Z(z) = D sin 



w z 
n 



, nTT c 

where w = — 

c n d 



Thus , 



the equation for Z(z), which gives the dependence of the acoustic 
variables on the layer depth, determines the eigenfrequencies w . 
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These are the cut-off frequencies of the various radially-symmetric 
normal modes of the layer. 

Use of the Laplace transform method to solve Eq. 2.2 with 
the boundary condition at the source written as 

u r U.O * llti'li*'), 

which implies excitation of only one mode, casts Eq. 2.2 into 
the simpler differential equation 



f _L i- _ / + 77 / \ r\ 



2.3 



where{6(r,s) is the Laplace transform of0(r,t). V(t) Z(z) is 

the radial velocity at the source. The relevant boundary conditions 

become 

4II ,VM 

Ji c I r - 0 . 



and 



< 2 $ -> O ocs r — *► 00 

The change of variable u = — 



c \ n 



2 2 

w + s casts Eq. 2.3 into the 



form 



1 - (?) - C> 



which has solution 



2.4 



(15) 



(jb - + ^ U> ) 



where I and K are the modified Bessel functions. Since I (u)-^OO 
00 o 

and K o (u)-^0 as u-*oo, the solution to Eq. 2.4 which obeys the 



required boundary conditions is 

<Z> -- CCl.) Ko u) 



where 
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Under the Laplace transform, the acoustic pressure, radial component 
of particle velocity, and radial component of particle displacement 

I3 0 p nm o 



where the contour is indented into the right half-plane around the 
singularities on the imaginary axis. Exact solutions to Eq. 2.5 
do not appear to be known, but certain approximations may be used 
to obtain approximate solutions for the leading portions of the 
propagating transient. 

For large values of argument, use of standard approximations 
for the modified Bessel functions results in the following asymptotic 
form 





Formal inversion results in expressions like 

Uoo 




2.5 



-too 






I n. ^ 1? i +* 



, > _Ttf / 



2.6 
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and 



cj 


and T I 


~~z~~ ) 




01 




JL a.\ 

6 ~) 






01 = 


+ 1 


<b Ox 

i~ 8 ru 


’ It 8 1 




11 = 




- 

a- / 






11 = 


AL. 

\ x *6 


I <6 

lie n- 


IS* 

a <3 < 


< A- 



The symbol T represents the time of flight over the shortest 
path between the surface of the source and the field point. The 
restriction on s necessary to satisfy the approximation which led 
to Eq. 2.6 is 

| 7 H ' i • 

Imaginary s greater than w^ in magnitude requires 



T- 5 Uk. 



where 



U)c*- 



- LOrv V 1 + ( O-uJa) 



is clearly a function of the source-radius and the cut-off frequency 
for the mode of interest. For a waveguide with pressure-release 



w a 
n 



boundaries which is excited in the lowest mode, the quantity 

Ta 

has the value which demonstrates the importance of the ratio 
of source-radius to layer-depth. As a result of the approximation 

which led to Eq. 2.6, Eq . 2.5 can be broken into two parts: 

Ur 2TT . . 

i(t1 *■ ‘ J 1 r L 

f “rr r C t 

* 3 v 

- iw^. ' 
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Each portion contributes significantly to u only when the exponent 
of each is stationary. For large s (which corresponds to the 
leading portion of the signal), only the first integral is important 
and yields its major contribution when the exponent has a vanishing 



derivative 

ci*- 



& O 



- 



(*.i - Ter) “ 0 - £ - ^ Ac 

Thus, the requirement that \bo/i \ ^ sets an upper limit 

to- -- T HTTiSf 

on the time within which the first integral yields the major contribution. 

For a step input, for which V(s) = the radial component of the 
particle displacement becomes 



£c_ 


1 


j 








K* ( T' 





Application of the asymptotic approximations of Eq. 2.6 leads to 

-Ter 



£ 






Inversion 



(5) 



leads to 



-\ 



where arg 



r (z\ n _ 4 

+ All ( ivi) 2 

= K n r? 



T . 



For a step input, the solution for the acoustic pressure becomes 

* • rz< { t / \ n f 



rr\r ^ 



3^-ia^v 6c, , 2 ' 8 
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The quantity t* = t - T appearing in Eqs. 2.7 and 2.8 represents 
the time which has elapsed since the propagating transient signal 
traveling with velocity c first arrived at the position r. This 



quantity will be denoted as the ‘’delayed time". 



For the leading portion of the propagating transient, one 
can apply a standard trigonometric asymptotic approximation^”*^ 

) 2 - 9 

to the functions in Eqs. 2.7 and 2.8. Application of this approx- 
imation to the sums of Eqs. 2.7 and 2.8 leads to expressions of 



the types 



J2- ivw a 



rx> r 









e cost 



r\ 






Application of Eqs. 2.9 and 2.10 to the solution for particle 
.acement casts it in the form 



2.11 



aw 

where 2TTN = y and N = %w T — — . 

y a n c 

The quantity N is the number of cycles of the transient generated 

in time t, and N corresponds to the cycle where the approximations 
a 

(and thus the theory) fail to apply. 

By combining the sine and cosine terms of Eq. 2.11, we find 
the expression for the radial component of particle displacement 
to be 



20 



2.12 



where C 



An -ti- 
tan ^ ~ “ ttiz. 



and 



1- Bm(N/W) Z ' 

Application of Eqs, 2,9 and 2.10 to the solution for acoustic 
pressure casts it in the form 



r-(¥in? mcKtoH-5*?) 



2.13 



where 



A(t) = | f 1** Bo* (nI/ n^ 1 ] 1 +■ L 



and 



tan 



va _ A&i (N / 

1- Bo\ (H/Ni*) 



V . . (5,14) 

Previous wort has shown that for the same input as above 

the exact solution for propagation of a transient in a duct becomes 
the leading term in the solution developed for the isovelocity layer, 
It should be noted that these duct terms can be retrieved from the 
trigonometric approximations by letting the amplitude factors go 
to unity and the phase angles go to zero. 

The solutions, using Bessel functions or asymptotic approxi- 
mations to the Bessel functions, for the transient response of an 
isovelocity layer have been developed with the restriction 

T - t - or Jsl/N*. — 1* 

Computer programs have been written to evaluate the predicted 
waveforms for acoustic pressure and particle displacement, using 
either the Bessel function or trigonometric approximations (see 
Appendices A and B, respectively). These programs give printed 
and graphical outputs for the duct and layer terms. In addition, 
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these programs predict the axis-crossings and amplitude-modulation 
in the layer. Figures 2.1 and 2.2 illustrate the acoustic pressure 
and particle displacement, respectively, predicted by the two approxi- 
mations. Tables 2.1 and 2.2 give the calculated values of the predicted 
axis-crossings and amplitude-modulation for acoustic pressure and 
particle displacement. It should be noted that the first portion 
of the first half- cycle of the waveform predicted by the trig 
computer program is not accurate. To avoid a mathematical ambiguity 
at the leading portion of the waveform, its value has been set equal 
to zero. For comparison purposes, the predicted amplitude peaks 
were normalized to the value of the highest peak in the waveform. 

As can be seen from a study of Tables 2.1 and 2.2, the difference 
between the two predictions is on the order of one per cent which 
is smaller than the experimental accuracy encountered. Therefore, 
the predictions from the computer program using trigonometric 
approximations will be used for all comparisons with experimental 
data . 

The approximate solutions for the propagation of a transient 
in an isovelocity layer have been developed, and computer programs 
to evaluate these solutions and predict the axis-crossings and 
amplitude-modulation in the layer have been written. 

Transfer Characteristics of Mylar Transducers 

To have some measure of confidence in the received waveform, 
it was decided to investigate the transfer characteristics, under 
steady-state conditions, for Mylar transducers. This investigation 
will cover the radiation impedance at the surface of the transmitter, 
r = a, and the radiation impedance of a receiver. 
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p, arbitrary units p, arbitrary units 

-CC3 OCG 003 -0C3 OGC 033 





Acoustic Pressure p vs. Delayed Time t' as Predicted by Computer 
Programs for the Conditions r=91.U cm, a=U.OO cm and d-8.00 cm. 

Figure 2.1 
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Table 2.1 



AXIS CROSSINGS 
USING BES 
TIME IN USEC 


AXIS CROSSINGS 
USING TRIG 
TIME IN USEC 


1.45 


1.39 


7.53 


7.46 


18.32 


18.26 


33.58 


33.51 


52.99 


52.91 


76.17 


76.08 


102.72 


102.64 


V 

132.33 


132.21 


164.57 


164.42 


199.13 


198.96 


235.72 


235.51 


274.07 


273.82 


313.96 


313.65 


355.18 


354.81 



AMPLITUDES OF 
MAXIMA, BES 
NORMALIZED 


AMPLITUDES OF 
MAXIMA, TRIG 
NORMALIZED 


1.000 


1.000 


0.632 


0.647 


0.470 


0.478 


0.390 


0.396 


0.339 


0.344 


0.303 


0.308 


0.275 


0.280 


0.253 


0.258 


0.235 


0.239 


0.219 


0.224 


0.205 


0.210 


0.193 


0.198 


0.182 


0.186 


0.171 


0.177 



Comparison of computed axis-crossings and amplitudes of maxima for 
acoustic pressure using Bessel functions and trigonometric approxi- 
mations to Bessel functions for the conditions r = 91.4 cm, a = 4.00 cm, 
and d = 8.00 cm. 
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Radial Component of the Particle Displacement vs. Delayed Time t' 
Predicted by Computer Programs for the Conditions r*91.8 cm, a=6.00 cm, 
and d«8.00 cm. 



Figure 2.2 
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Table 2.2 



AXIS CROSSINGS 
USING BES 
TIME IN USEC 

3.75 

12.48 

25.94 

43.82 

65.75 

91.33 

120.16 

151.85 

186.05 

222.44 

260.73 

300.68 

342.07 



AXIS CROSSINGS 
USING TRIG 
TIME IN USEC 

3.92 

12.60 

25.99 

43.79 

65.64 

91.13 

119.89 

151.51 

185.65 

221.99 

260.24 

300.15 

341.50 



AMPLITUDES OF 
MAXIMA, BES 
NORMALIZED 

0.474 

0.701 

0.844 

0.933 

0.982 

1.000 

0.995 

0.974 

0.941 

0.901 

0.857 

0.811 

0.764 



AMPLITUDES OF 
MAXIMA, TRIG 
NORMALIZED 

0.475 

0.698 

0.841 

0.931 

0.981 

1.000 

0.996 

0.975 

0.943 

0.903 

0.860 

0.814 

0.768 



Comparison of computed axis-crossings and amplitudes of maxima 
for the acoustic particle displacement using Bessel functions and 
trigonometric approximations to Bessel functions for the conditions 
r = 91.8 cm, a = 6.00 cm, and d = 8.00 cm. 
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In a steady-state condition, the acoustic potential for a 



cylindrical source in a waveguide can be expressed as 
The time independent wave equation becomes 








0 < 



2.14 



The separation of variables suggested in Eq. 2.14 leads to the 

value of the separation constant, k = — , where the phase velocity 

P 

c is defined as 

P 

c ? ji- (“~F 2,15 



Thus, the wave equation for the radial component of the acoustic 



potential becomes 



£1 , x v I = o 

Jr 1 r Ar 



with the general solution 



“ A.”^o ^ • 2.16 



For an outward traveling wave, Eq. 2.16 becomes 

R (r) - Afjo d] = AlWkVr). 



2.17 



The acoustic potential, acoustic pressure, and the radial component 
of the particle velocity become 
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Using Eq. 2.15, the specific acoustic impedance £ becomes 

Z. = ^ i(>c 



where arg = 



wr 



fPTW- 



2.18 



Then, the radiation impedance Z is given by 

Mo* 1 '* 






ir-Ki. - i P c P A H ° ^ 

r - ' T I4,<“ 

where A is the total radiation area of the transducer. Thus, the 
radiation impedance at the surface of the transducer, r = a, is 
given by 

2 _ j. P cA He 0 ^ il~ (*&? ) 

r = fnlf 7WF) ‘ 2 , 9 



$ 



The behavior of Eq. 2.19 as a function of frequency will be 
investigated. An estimate of the maximum frequency encountered 
in a shallow-water layer can be derived by recalling that the 
highest frequencies in the transient occur in the first half-cycle. 
The following relation concerning the time for the maximum frequency 
can be derived: 



W^tw-T* - Z.T1 N -TT or 

- T 2, + 1 



The relation for the instantaneous frequency 
can then be manipulated into the desired result 









-\j i t . 
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For r-a = 100 cm and d = 8 cm, 

w = 650 khz. 
max 

The cut-off frequency for the first mode and the above approximations 
is 



o 

w =59 khz. 
n 

Thus, the range of frequencies expected to be encountered are 
v/^ 4 w i Wmax* The expected frequency must always be greater than 
the lower limit, for the assumption of a steady state condition no 
longer holds at w = w^. 

Since the Hankel functions are complex numbers, they can be 
expressed as a modulus and a phase angle. Thus, Eq. 2.19 can be 
rewritten as 



A 




e. 



where Mod is the modulus and 0 is the phase angle. Figure 2.3 
shows the dependence of the phase angle and modulus on the argument 
of the Hankel functions. For small values of the arguments, the 
radiation impedance is strongly dependent on the frequency. How- 
ever, when the value of the argument is about IT which corresponds 
to a frequency within 0.5 per cent of the cut-off frequency, the 
radiation impedance can be approximated by A^c . Thus, for the 
frequency range of interest, the radiation impedance is frequency 
independent . 

The open-circuit mechanical impedance of a Mylar transducer 
is given by = r + i(wm - s/w) where r is the transducer mechanical 
resistance and s is the transducer mechanical stiffness. For the 
frequency range of interest, the mechanical resistance is considered 
to be negligible in comparison with the reactive part of Z^. Under 
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3 
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the Argument of the liarikel Functions 



/ 1 c \ 

the assumption that the mechanical stiffness is given by 
s = P A/jj) where P q is the atmospheric pressure and ^ is the 
effective depth of the air cushion behind the Mylar diaphragm, the 

equation for Z can be rewritten as 

m -- 

^ L A L w f wt ^- 

where t is the thickness of the Mylar, With the nominal values 
pmylar = U395 m/ ™' 

-3 

t = 2.54 x 10 cm, 

P q = 1.013 x 10^ dyne/cm^, and 
j(= 2,00 x 10" 3 cm, 



the values of for the frequency range of interest can be 

calculated. Figure 2.4 shows the dependence of the magnitude of 
Z 

— vs o frequency. 

The velocity of the diaphragm of a Mylar transducer is 
given by 



i iO v 

-u) z ECo 



2.20 



where is the blocked diaphragm electrical capacitance, v is 
the signal voltage across the transducer, 0 is the turns ratio, 
and R is the internal resistance of the power amplifier. The 
turns ratio 0 is defined as 

0 -- 

where V q is the external polarizing voltage and x q is the position 
of the Mylar film when the polarizing voltage is applied. A 
comparison of the magnitudes of the terms in the denominator of 
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Modulus of Motional Impedance Z m of the Mylar Transducer used as a Receiver vs. Frequency. 

Figure 2.k 
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Phase Angle and Mod vs. Frequency for i v jylar Transducer used as a Receiver. 



is proportional to the particle displacement. These conclusions 
only hold for the leading portions of the received transient where 
can be replaced by A|>c. 

The effect of the transfer characteristics of Mylar transducers 
has been found to be dependent on the frequency and the radius of 
the transducer. Given a transient signal which is the mirror of 
the input voltage, the following is a summary of the effects that 
would be observed as a result of the radiation impedance: 

(a) since the radii of the family of transmitters were large, 
the response is essentially flat throughout the region of interest; 
thus, no effects should be observed. 

(b) because of the small radius of the receiver, the effect 
becomes more pronounced than for the transmitters. One should, 
observe a reduction in the size of the peaks of the leading portion 
of the transient and should observe an enhancement of the peaks as 
the local apparent frequency approaches the cut-off frequency. A 
similar phase angle effect should be observed. 
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3. EXPERIMENTAL APPARATUS 



The experimental apparatus in this research project consisted 
of a tank, a family of cylindrical Mylar transducers, a barium 
titanate array receiver, and the associated electronics equipment. 
The laboratory setup is shown in Fig. 3.1. 

Tank 

The water layer was contained in a tank constructed of 5/8 
inch marine plywood. The dimensions of the tank were 6 ft x 6 ft x 
1 ft. The plywood was covered with Varathene to provide water- 
proofing and lined with 1/4 inch polyethylene foam which was assumed 
to provide the desired pressure-release bottom boundary of the 
layer. The polyethylene, not being waterproof, was coated with 
liquid Neoprene paint to prevent deterioration of the pressure- 
release surface. The bottom of the water layer was maintained 
level by means of adjustable feet installed in the legs of the 
supporting platform. Since the walls of the tank were also lined 
with polyethylene, the tank could be excited as a pressure-release 
cavity, which provided one method of determining the depth of the 
layer . 

Transducers 

Previous work^’^) had indicated that Mylar dielectric 
transducers were well- suited for use in layers. 

A family of eight Mylar transducers and one barium titanate 
linear array were constructed for use in this investigation. The 
radii of the Mylar transducers were 10 cm, 6 cm, 4 cm, 2 cm, 1 cm, 
0.5 cm, and two at 0.25 cm. Each had an active radiating face 
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8 cm in height which was the maximum depth of the layer. These 

various transducer dimensions were chosen to provide a large range 

of values for the ratio N/N . The barium titanate array was 8 cm 

a 

high and 1/8 inch in diameter. 

The Mylar transducers were constructed using techniques 
/ 1 £ \ 

described in , modified as follows: 

(a) The backplates of the larger transducers were constructed 
from \ inch thick aluminum cylinders. The backplates for the 
smallest transducers were fabricated from \ inch diameter copper 
tubing. 

(b) The cylinders were capped at the top and bottom by 
lucite plates, 

(c) The Mylar film was cut to the desired dimensions and 
attached to the aluminum cylinders with thin strips of double-stick 
tape, one piece around the top, one around the bottom, and one 
piece vertically from bottom to top. 

(d) The inner electrical connection was made from the center 
pin of a BNC connector mounted in the center of the top lucite 
plate to the inside of the backplate. The outer conducting path 
from the metallized side of the Mylar film to the BNC was completed 
by covering the top lucite plug with aluminum foil which was lapped 
over the Mylar film* To ensure good electrical continuity, silver 
paint was applied to the aluminum-Mylar interface. To decrease the 
possibility of electrical pickup, the bottom lucite plate was also 
covered with aluminum foil. 

(e) When the silver paint had dried, the entire transducer 
was covered with liquid neoprene to provide the necessary water- 
proofing. 
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Tank with transmitter and receiver. Electronics equipment in 
the background. 

Laboratory Setup 
Ei gnre 3.1 




Transmitters and Receivers. 

Ei gure 3 . 2 
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The pressure-sensitive barium titanate was constructed using 
eight 1/8 inch by 1/8 inch cylinders mounted on several strands 
of twisted wire, which provided the inner electrical connection. 

To provide the damping necessary to reduce the effect of longitudinal 
resonances, the spaces between the cylinders were filled with strips 
of rubber. The outer electrical connection was provided by a wire 
which was coiled around each cylinder. This probe proved to be 
somewhat insensitive and was used only for comparison purposes. 

When the electrical load on the Mylar receiver (see Fig. 2.5) 
was terminated by the 1 megohm impedance of the voltage amplifier, 
the receiver was sensitive to particle displacement. Therefore, 
the Mylar receivers were used to measure the particle displacement, 
and the barium titanate array was used to measure the acoustic 
pressure. Transmitters and receivers are shown in Fig. 3.2. 

Electronic Equipment 

Figure 3.3 shows the block diagram of the equipment; Table 
3.1 identifies the components in the figure. The unit-pulse generator 
was used for the step input. 

The received signal was sent from the receiver through the load 
and bias network, amplifier and filters to trace "B" of the dual- 
beam oscilloscope. The output of the power amplifier could be 
displayed on trace "A". Measurements of the received waveform axis 
crossings and amplitudes were made directly from the oscilloscope 
display, or photographed for later analysis. The electronics were 
triggered by the time-mark generator which also triggered the time- 
base of the oscilloscope. The intensity of the "B" trace was 
modulated by a 10 usee marker available from the time-mark generator. 
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Table 3.1 





Electronics Equipment 


UPG 


Unit pulse generator, General Radio type 1217-B 


PWR AMP 


Power amplifier, Hewlett Packard 467-A 


TMG 


Time Mark Generator, Tektronic type 180A 


AMP 1 


Voltage amplifier, Tektronic type 1121A 


SKL FIL 


Electronic band pass filter, SKL model 202-D 


RC FILl 


Homemade variable low-pass RC filter with 400 khz roll-over 


RC FIL2 


Homemade high-pass RC filter with 2 khz roll-over 


SCOPE 


Dual beam oscilloscope, Fairchild model 777 


LBN 


Load and bias network shown in Fig. 3.4. 
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Source and receiver were carefully positioned so that the distance 
between adjacent sides, r-a, was always the same* This insured 
that the trace modulation would be properly positioned with 
respect to the inception of the received transient signal. The 
load and bias network shown in Fig. 3.4 served to supply the D.C. 
voltage bias needed by the Mylar transducers without allowing 
this voltage to appear elsewhere in the circuit and to provide 
the proper electrical load for the receiver. The 22 ohm resistor 
in the transmitter circuit was needed to suppress parasitic 
oscillations caused by the transient excitation of the capacitance 
of the transmitter and the residual inductance of the power 
amplifier. 

The rise time of the system was experimentally determined 
to be about 0.5 usee. 
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4, EXPERIMENTAL INVESTIGATION 



The purposes of the experimental investigation were to (a) observe 
and analyze the received acoustic pressure and radial particle displace- 
ment waveforms in the layer, (b) compare these waveforms with the 
theoretical predicted waveforms produced by the computer programs, and 

(c) vary the value of the ratio N/N over a large range to investigate 

a 

the validity of the assumed asymptotic approximations previously dis- 
cussed, Four basic parameters characterized each experiment: the 

range r, the transmitter radius a, the layer depth d, and the speed 
of sound, c. 

Experimental Procedures 

Before each run, it was decided which parameter would be varied. 
Generally, the time of flight T and the depth of the layer d were 
kept constant while the source radius a was varied throughout its 
range. The 0.25 cm radius Mylar transducer was used to record the 
received particle displacement, and the barium titanate linear array 
was used for the acoustic pressure receiver. The temperature of the 
water was measured to determine the speed of sound in the layer . 

The depth was determined by averaging several physical measurements 
taken over the ensonified area. Another method of determining the 
depth was to excite the tank as a pressure-release cavity. The 
cut-off frequency could then be determined by observing the received 
signal on an oscilloscope and varying the frequency until the signal 

disappeared. Once the cut-off frequency was known, the depth could 

TTc 

be calculated, using the formula w^ = — * Because of the irregularities 

observed in the depth, this method proved to be no more accurate than the 
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actual measurement using a ruler and was consequently abandoned. 
The intensity of the trace of the oscilloscope was modulated by 
a 10 usee marker available from the time-nark generator. The 
source and receiver were carefully positioned so that the distance 
between adjacent sides, r-a, was always the same and such that 
the reception of the received transient and one of the time marks 
always coincided. This insured that when changing transmitters, 
identical positioning and time of flight could be achieved. 

The received waveforms were photographed from the oscil- 
loscope display. For the step input, several photographs were 
taken of each waveform. An expanded time scale was used when 
photographing the leading portion of the waveform which contained 
the high frequencies. For a typical experimental run, three 
photographs of each waveform were made utilizing, for example, 

20 usec/cm time scale. By varying the delay time vernier, the 
entire waveform could ba displayed and photographed for later 
analysis. By careful measurement of the position of each axis 
crossing with the neighboring time marks, the relative phase of 
each axis crossing could be compared between the photographs taken 
for different source radii. The amplitudes of the peaks of the 
waveforms were read with the help of dividers and a scale. Sample 
photographs are shown in Fig. 4.1. It should be noted that the 
anomalous behavior of the received waveforms as evidenced by 
Fig. 4.1 departs drastically from the smoothly- vary ing computer 
predictions of Figs. 2.1 and 2.2. This anomaly will be discussed 
in more detail later. 
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Step Input 
in Velocity 

r=63.0 cm 

a=6.00 cm 

d=3.00 cm 

temp=22°o 



l, a) Received Acoustic Pressure 




Step Input 
in Velocity 

r= 91 .l 4 cm 

a=6.00 cm 

d=8.00 cm 

tenp=l6.5 C 



Observed Waveforms 
'•’i.'”. -e . 1 . 
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For each experimental run, a corresponding computer run was 

made using the same values for r, a, c, and d. The computer 

programs were written so that their output for each run included 

graphs of the predicted waveforms, times of axis crossings for 

the layer solution t* and for the duct solution t' , and amplitudes 

m on 

of the peaks and the time of occurrence of th^se peaks for the 

layer solution. In this manner, th 3 theoretical values of the 

times of axis crossings and of the amplitude peaks could be 

compared with those values experimentally determined. 

The time of the axis crossings were compared in the following 

manner: Previous work^^ suggested that the time of each axis 

crossing should be the same as the theoretically predicted t* 

rr. 

within the validity of the approximations used. As was shown in 

section 2, the theory was valid only as long as N/N £ 1. There- 

a 

fore, the convenient parameters for analyses were the time of the 
axis crossings and the ratio N/N^ . In ord^r to reduce the size of 

cl 

the graphs, the time coordinate was reduced by the amount t^ predicted 

by the computer.. The time coordinate signified by t 1 was plotted 

against the ratio N/N . Figures 4.2 to 4.4 are representative 

a 

examples of these graphs. The solid line represents the values of 

t* - t', and the circled dots represent the values of t * - t 1 . 

om m exp m 

Were the theory to have been completely correct and the experiment 

accurate enough, the data would have fallen on the X-axis (which 

represents t* - t'). As can be seen from Figs. 4.2 to 4.4, the 
mm 

data were not in agreement with the theory. The effect of the 
phase angle demonstrated in Fig. 2.6 was taken into account and 
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Experimental and Computed t 1 vs. N/N for Acoustic Pressure Waveform 
Generated by a Step Input in Velocity. 



Figure it. 2 
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Experimental and Computed t* vs. N/Nn for Particle Displacement 
Waveform Generated by a Step Input in Velocity for a = 1.0 cm. 

Figure U.3 
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Experimental and Computed t^ vs. N/K a for Particle Displacement 
lave form Generated by a Step Input in Velocity for a = 6.0 cm. 



Figure 1+ 
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plotted on Figs. 4.3 and 4.4. As can be seen, this did not 
improve the agreement between the theory and the data. 



The peak- amplitude data were also analyzed. The experimental 
amplitude peaks were first normalized to the value of the highest 
peak in the waveform. The predicted peaks were normalized in the 
course of the computer programs. The ratio of these normalized 
quantities was then plotted as a function of the delayed time of 
the predicted amplitude peaks for each run. Exact agreement in 
the shape of the transient would result in a line of zero slope 
passing through the point (1,0). ‘The amplitude effect discussed 
in Section 2 for a Mylar transducer used as a receiver was also 
plotted. Figures 4.5 to 4.7 show plots for various values of r, 
a, and d. Once again, this effect did not improve the agreement 
between the theory and the data. 

Error Analysis 

The rise time of the experimental system was determined to be 
about 0.5 usee. Therefore, instead of being a perfect s tep- function 
input, the velocity at the source was initially described by 



should introduce distortion into the received signal for at most 
3.0 usee of delayed time t*. Since this t ime is less than that 
for the first maxima, the effect of the rise time on the received 




where o( = 2 x 10^ sec ^ is the reciprocal of the rise time. 
Previous work^^ has indicated that the rise time of the system 



signal can be ignored in this investigation. 
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a Step Input in Velocity vs. Delayed Time for a = 1.0 cm. 



Experimental Results 
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Figure 4.7 



It has been shown that the velocity of the Mylar face decays 
to 88 per cent of its maximum value after a time of 360 usec^\ 

This effect is not Large enough to account for the tailing off 
of the amplitude peaks of the received waveforms seen in Figs. 

4.5 to 4.7. 

Uncertainties in r, a, c, and d were the possible sources 
of error in this investigation. The effective radius of each 
transmitter was measured to the nearest 0.5 mm, which could have 
introduced an uncertainty of + 1.0 per nt in the value of a 
for the 6.0 cm radius source and +5.0 per cent in the value of a 
for the 1.0 cm radius source. Measurement of the temperature of 
the water to the nearest half-degree restricted the systematic 
error in c to within + 0.5 per cent. Determination of the time- 
of- flight T of the received signal form the source to the receiver 
produced a possible systematic error of + 0.5 per cent, which in 
turn resulted in an error in r of + 0.5 per cent. 

The largest experimental error was encountered in determining 
the depth of the layer. Observed irregularities of + 1.0 mm in 
the bottom surface of the tank could introduce an error in the 
local effective depth of the layer. Since the averaged layer depth 
used as an input to the computer programs was only accurate to + 1.0 
mm, the systematic error introduced in w was on the order of + 3.0 
per cent for a depth of 8.0 cm. 

The error in determining the times of the axis crossings from 
the data was effected by the accuracy of the oscilloscope time-base 
which was calibrated to within 0.3 per cent. With the trace of 
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the oscilloscope modulated with the 10 usee marker as an additional 

time reference, this error was reduced even more. Therefore, a 

random error of + 1.0 per cent could be expected when one takes 

into account the accuracy of reading the oscilloscope scale. 

The errors in the computed values of t‘ were influenced 

m 

predominately by the error in w , since it appears in the argument 

r 2 ' T 

of the Bessel functions (w J t - T ) from which the t' are 

n V 7 m 

computed. It has been shown^ that the argument of the Bessel 
functions could be rewritten as 






where C are a set of constants which yield zero values for the 
m 

amplitude of the particle-displacement solution. For the leading 
portion of the signal, we have 



U) 



* {rr 



c 






and 



c 1 

V- KV\ 



t'„ ~ ut 2.T 

Therefore, the +3.0 per cent error in w could lead to a + 6.0 
— n — 

per cent error in the computed values of t\ 

To reduce the scatter of the data, a one-parameter fit was 
tried. It was decided to determine an effective cut-off frequency 
and use this as an input to the computer programs instead of the 
layer depth. This was accomplished by determining the cut-off 
frequency for the values of t anc * t *^ ie arguments of the trig 
function approximations developed in Section 2. At an axis 
crossing, the argument was equal to an integer multiple of IT 
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A computer program was written to evaluate the cut-off frequency 
for each axis crossing. Then, these cut-off frequencies were 
averaged to give an "effective cut-off frequency" to be used in 
the computer programs to recompute the time for each axis crossing. 
Figure 4.8 shows the effect of this one- parameter fit. As can be 
seen, it did not improve the scattering significantly and was not 
used . 

Error flags on Figs. 4.2 to 4.8 indicate the relative magnitudes 
of errors in the experimental and computed values of the axis crossings 
and peak amplitudes. However, these error flags are not large enough 
to explain the scatter of the data or the anomaly which was present 
(see Fig. 4.1). The theory could not account for this anomaly so 
various experiments were conducted to investigate its character and 
source of origin. 

Firstly, it was suggested that irregularities in the bottom of 
the tank might be a possible cause of the anomaly. Therefore, the 
relative positions of the transmitter and the receiver in the tank 
were changed. The values of T and d were always kept constant. Four 
contrasting positions were chosen; one near the left hand side of 
the tank, one in the center of the tank, one adjacent to a seam in 
the bottom near the right hand side, and one where the transmitter 
and receiver straddled the seam. A.s can be seen in Fig, 4.9, within 
the experimental accuracy of the investigation, positioning on the 
bottom had no effect on the received waveform. 

Secondly, it was decided to monitor the current drawn by the 
transmitter when loaded by the water medium and compare this with 
the current predicted assuming pc loading. The current was 



58 




59 



Experimental and Computed t^ vs. M/N a for Particle Displacement V.aveform Generated by 
a Step Input in Velocity using a One Parameter Fit. 

Figure U . 8 
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experimentally determined by measuring the current in the two legs 

of the circuit shown in Fig. 4,10. The current present in each 

leg was amplified by a differential amplifier and then the difference 

displayed on the oscilloscope. Photographs of the received signal 

could then be analyzed to determine the current drawn as a result 

of acoustic radiation into the medium. By analyzing the equivalent 

a6 x 

circuit of a Mylar transducer used as a transmitter ^ , as shown 

in Fig. 4.11, the predicted value of the current drawn as the result 
of the assumed pc loading can be determined. For this assumed 
loading, the current drawn by the mechanical leg of the source 
should be flat. Typical values for the 6 cm radius source were 
near 6 uamp. The measured values of the current had rapid oscil- 
lations of nearly 1 uamp peak-to-peak amplitude, which indicated a 
significant fluctuation from the flat response predicted. This 
casts doubts on the validity of the assumption of pc loading used 
in the developments of Section 2. 

Lastly, it was suggested that the polyethylene foam was not a 
perfect pressure-release surface and that transmission through it 
into the wood bottom was the cause of the anomaly. In order to 
investigate this possibility, two additional layers of polyethylene 
were glued to the bottom of half the tank. Then runs were made on 
the s ingle- layered bottom and then on the three- layered bottom. 

The quantities r, a, and d were kept constant for the corresponding 
runs. Figures 4.12 to 4.14 indicate the predicted waveform, the 
experimentally determined waveform on the three- layered bottom, and 
the experimentally determined waveform on the single- layered bottom 
for three values of a: 10, 6, and 4 cm. The vertical and horizontal 
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Circuit used to Monitor Current Drawn by a Mylar Transducer as a 
Result of Acoustic Radiation into a Water Medium. 

Figure U.10 




6 2 

c = 1.5 x 10 kg/sec-m v(t) = 20 volt step input 
Equivalent Circuit of a >tylar Transducer used as a Transmitter. 

Figure lull 
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Step Input 
in Velocity 

r=70.66 cm 
a=10.00 cm 
d=8.00 cm 



(b) Received particle displacement waveform on three -layered 
bottom 




Step Input 
in Velocity 

r=70.66 cm 
a=10.00 cm 
d=8.00 cm 



(c) Received particle displacement waveform on single-layered 
bottom 

Bottom Layer Investigation for Three-Layered Bottom and Single-Layered 
Bottom using 10.00 cm Source. 



Figure U.12 
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arbitrary units 




(a) Predicted Waveform 



Step Input 
in Velocity 

r-66.66 cm 
a-6.00 cm 
d=8.00 cm 




. (b) Received particle displacement waveform on three-layered 
bottom 




Step Input 
in Velocity 

r-66.66 cm 
a=6.00 cm 
d=8.00 cm 



(c) Received particle displacement waveform on single-layered 
bottom 

Bottom Layer Investigation for Three -Layered Bottom and Single-Layered 
Bottom using 6.00 cm Source. 



Figure U.13 
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(b) Received particle displacement xjaveform on three-layered 
bottom 




Step Input 
in Velocity 

r=6ii.66 cm 
a=U.OO cm 
d=8.00 cm 



(c) Received particle displacement waveform on single-layered 
bottom 

Bottom Layer Investigation for Three-Layered Bottom and Single-Layered 
Bottom using U.00 cm Source. 



Figure U.lU 
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scale settings on the oscilloscope were kept the same for each 
set of runs* As can be seen, there is some evidence of a bottom 
effect, but there is no clear correlation between the three runs. 
However, this evidence is sufficient to doubt the credibility of 
the assumption of the perfect pressure-release characteristics of 
the 0*25 inch polyethylene foam. 
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5. RESULTS AND CONCLUSIONS 



The investigation has shown that for N/N £ 1 the predicted 
and experimentally observed shifts in axis crossings should agree 
to within + 6 per cent if one takes into account just the known 
systematic errors. The correlation between the predicted and 
observed amplitudes should be of the same accuracy. However, 
the difference between the theory and data is greater than that 
expected , 

The theory as developed in Section 2 could not account for 
the anomaly which was present in all runs, but the qualitative 
agreement between the observed and predicted waveforms and the 
general agreement with the solutions developed in other work^ 
lead one to believe that the theory presented in Section 2 is correct. 

It is concluded that the experimental setup used is not accurate 
enough to give results needed to confirm the theory. 

It is recommended that better electronic equipment and transducers 
should be investigated. Further studies should be conducted to check 
the validity of the assumed pressure-release characteristics of the 
polyethylene foam. More investigations are needed to try to determine 
the cause of the observed anomaly. 
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AMPXP=FXPER I MENTALLY DETFF MI NED VALUES PF PEAK 

amplitudes <* on -) in pressure. 
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N2=NIJMHFR of EXPERIMENTAL maxima (♦ OR -) IN 
PRESSURE. 

AMPXP=EXPFR IMENTAU.Y DETFPmimeo VALUFS OF PFAK 

amplitudes ( + or -» in prfssure. 
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